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ABSTRACT 

We present an 8a detection of cosmic magnification measured by the variation of quasar density 
due to gravitational lensing by foreground large scale structure. To make this measurement we used 
3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing 
~ 200, 000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation 
function exhibits the amplitude, angular dependence and change in sign as a function of the slope of 
the observed quasar number counts that is expected from magnification bias due to weak gravitational 
lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, 
seeing variations and errors in the photometric redshifts) are well controlled and do not significantly 
affect the lensing signal. By weighting the quasars with the number count slope, we combine the 
cross-correlation of quasars for our full magnitude range and detect the lensing signal at > 4cr in all 
five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60/i _1 kpc to 
10 h^ 1 Mpc and are in good agreement with theoretical predictions based on the WMAP concordance 
cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide 
useful constraints on the galaxy-mass power spectrum. 

Subject headings: cosmology - large-scale structures - gravitational lensing: magnification - quasars 
- galaxies 



1. INTRODUCTION 

We expect the large-scale structure seen in the low 
rcdshift Universe to gravitationally lens background 
sources, such as high redshift galaxies and quasars. This 
lensing effect causes both a magnification and a distor- 
tion of these distant sources. The systematic distortion 
of faint background galaxies by gravitational lensing, 
the cosmic shear, has now been measured by several 
groups in the past few years llVan Waerbeke et al.ll2000t 
IBacon et ahl 120001 IRhodes et all 120011 IHoekstra, et al.l 
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Brown et all 120031 iMassev et all I2004|l and has been 
found to be in remarkable agreement with theoretical 
predictions based on the Cold Dark Matter model. 
It has also provided new constraints on cosmological 
parameters, especially on o~%,£l m and the shape of 
the dark matter power spectrum (for a review, see 
Van Waerbeke fc Mellierll2003t iRefregienEM iHoekstral 
2003). In addition to shear-shear correlations, the 
cross-correlation of foreground galaxies with back- 
ground shear, known a s galaxy-galaxy lensing has 
also b een measured llBrainerd. Blandforo^^Sjn^ii 
19961 IdelF Antonio fc Tvsonl Il996t iGriffiths et alJ 
T99a iHudsonet al l 11993: iFischer et all I200C 



Wilson et all 120011 i McKay et all 120011: iSmith et all 



20011 IHoekstra et alJ 120031) . Recent measurements 
from the Sloan Digital Sky Survey (SDSS; lYork et alJ 
2000) have enabled accurate constraints on galaxy halo 
profiles or more generally the galaxy- mass correlation 
ijSeliak et al1l2lMlSheTdon et al.ll200l . 

In a similar way, the systematic magnification of back- 
ground sources near foreground matter over-densities, 
the cosmic magnification, can be measured and can pro- 
vide largely independent constraints on cosmological pa- 
rameters. Gravitational magnification in the weak limit 
has two effects: First, the flux received from distant 
sources is increased, resulting in a relatively deeper ap- 
parent magnitude limited survey. Second, the solid angle 
is stretched, diluting the surface density of source images 
on the sky. The net result of these competing effects 
is an induced cross-correlation between physically sepa- 
rated populations that depends on how the loss of sources 
due to dilution is balanced by the gain of sources due to 
flux magnification. Any type of background source can 
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be used to measure this effect (galaxies, quasars, super- 
novae, etc.), but in practice, previous investigations have 
used foreground galaxies and background quasars moti- 
vated by the large redshift range probed by quasars and 
general redshift segregation between these two popula- 
tions. Despite the apparent elegance of this solution, 
lensing-induced quasar-galaxy correlations have been a 
controversial subject for more than a decade. Numerous 
teams have attempted to measure this effect and have re- 
ported detections of changes in the density of background 
quasars in the vicinity of galaxies. However, as seen by 
reviewing the literature on this type of measurement, the 
results have been generally discrepant with each other as 
well as in disagreement with the expected signal from 
gravitational lensing. 

The fir st analysis of qua s ar-gal axy correlations was 
done by iSeldner fc Peebles! (|10T9f) using a sample of 
~ 400 quasars and galaxies from the Lick catalog that led 
to a 3.7a detection of an quasar excess on ~ 20' scales in 
the vicinity of galaxies. The first measurements aimed at 
detecting the expected lensing signal used radio-selected 
quasar samples; this method yielded numerically larger 
quasar samples as well as a st eeper number cou nt relation 
to enhance the lensing signal. iFuemannl (^990) correlated 
bright, radio- loud quasars at moderate and high redshifts 
with galaxies from t he Lick catalog and found an ex- 
cess on a 10' scale. iBartehnann fc Schneider! (jl993b) 
repeated the analysis with 56 z > 0.75 optically iden- 
tified quasars from the 1-Jansky catalog and confirmed 
the previous result. Similar excesses were also found by 
cross-corr elating the 1-Jansky quasar catalog with IRAS 
galaxies (Bartclma nn fc Schneider! 119941: iBartsch et al. 

t and dif f use X -ray emission ( Bartclmajmj^LaL 
iCooravl Il999[) . iRodrigues- Williams fc Hogari 
found a correlation between optically-selected 
quasars and Zwicky clusters, but with an ampli- 
tude that cannot be reproduced by lensin g of sim- 
ple mass models. iSeitz fc Schneider! l|1995f) revisited 
the previous 1-Jansky/IRAS analysis, finding agree- 
ment for the intermediate redshift quasars but fail- 
ing to detect any corre lation for the high redshift 
ones. IWu fe Hani l)1995f) repeated the cluster cross- 
correlation using Abell clusters and found no correla- 
tion with the 1-Jansky sources . Using wide-field R- 

band images, iNorman fc Impevl l)1999|) detected a cor- 
relat ion with 1-Jansky quasa rs on scales greater than 
10'. iWilliams fc Irw in (199l| and INorman fc William^ 
(2000) cross-co rrelated LBQS and 1 -Jansky quasars with 
APM galaxies ijMaddox et alJll990]) and claimed signifi- 
cant over-densities on angular scales of the order of one 
degree. 

Using optically selected sources, similarly mixed re- 
sults have been obtained: by cor relating UV-excess 
quasars and APM galaxies in clusters. iBovle et all l|1988f) 
found a 30% deficit of qu a sars o n scales of 4' around 
galaxies. iCroom fc Shanks! l)1999|) investigated the lens- 
ing explanat ion but found the amplitude of the sig nal to 
be too high. IRodrigues- Williams fc Hawkins! l) 19951) used 
variability-selected quasars and found a correlation with 
Zwicky clusters that they interpret as induced by lensing. 
Again, the amplitude of the correlation largely exceeded 
expected results. Associations between red galaxies from 
the APM catalo g and moderate-redshift quasa r s wer e 
investigated by iBenftez fc Martinez-Gonzalez! fl995). 



They reported a ~ 30% excess of quasars within 2' from 
the gal axies and a signal con sistent with zero on larger 
scales. iFerreras et al.l l)1997|) cross-correlated optically 
selected bright quasars and galaxies. The amplitude and 
the redshift dependence of their results was inconsistent 
with either the lensing or the dust explanation, and they 
suggested that the quasar selecti on process su f fered f rom 
incompleteness. More recently, iMvers et~aH (|2003^ in- 
vestigated correlations between galaxy groups and opti- 
cally selected quasars from the 2dF survey. They found 
a 3<7 anti-correlation within 10', whose amplitude im- 
plies that the velocity dispersion of galaxy groups is of 
the order of 1000 km s -1 , i.e. much higher than ex- 
pected. In addition, they have investigated the effects of 
extinc tion by dus t and found them to be negligible. Fi- 
nally, Gaztanaga (2003) measured the cross-correlation 
between photometric galaxie s and spectroscopic qu asars 
using only t he SPSS EDR llStoughton et all 12002) . In 
contrast to IMvers et al! (|2003ft thev found a positive 
cross-correlation of 20% on arcminute scales - two re- 
sults that might be complementary because they probe 
quasar samples of different apparent magnitude (see, e.g., 
IMvers et alll2005|) . although the amplitude of both mea- 
surements was far in excess of the expe cted lensing signal . 
Compared to our measurements, the ( laztaiiaga (2003) 
analysis was done on earlier reductions of EDR data and 
an incomplete sample of quasars, both of which may have 
affected the observed signal (see {§yjfor more details). In 
addition, the foreground and background samples used 
appear to have significant redshift overlap, which can 
lead to a much stronger non-lensing correlation. 

Clearly, the scatter in the existing observational re- 
sults is large, ranging from significant positive correla- 
tions to null and negative correlations, as well as a vari- 
ety of claimed scales for the different detections, from half 
an arcminute to one degree. In addition, quasar-galaxy 
correlations have been controversial for many years due 
to the large discrepancy between the claimed detections 
and early theoretical estimations llSchneider et alJll992t 
| Bartelmannll995HDolag fc Bartelmannll997tlSanz et alJ 
1997). Indeed, the observed amplitude of previous mag- 
nification results was typically an order of magnitude 
larger than predicted by theory. Initially, it was not clear 
whether this problem was due to the observations or the 
models, as the early formalism describing cosmic mag- 
nification used several simplifying assumptions: a lin- 
earized magnification and a constant bias for the galaxy- 
dark matter relation. Improvemen ts have recently been 
made on the theoretical side, i.e.. iMenard et al. (2003) 
went beyond the linearized magnification approximation 
by including non-linear corrections to the relation be- 
tween the magnification and density fluctuations, and 
then compared their re sults to numerical simulations. 
iGuimaraes et, all l)2001f ) introduced a scale-dependent 
galaxy bias obtained fro m measurements of galax y auto- 
correlation functions. Uain. Scranton fc Shethl l)2003f) 
modeled the complex behavi or of the galaxy bias using 
the halo-model approach, and Takada & Hamana ( 2003) 
added an estimation of the full non-linear magnifica- 
tion contrib ution by including the mag nification profile of 
NFW halos l|Navarro. Frenk. fc White!ll997l) in the halo- 
model formalism. Together, these works provide an ac- 
curate theoretical framework and show that earlier es- 
timations of the amplitude of the cosmic magnification 
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were underestimated by 20 to 30 %. However, this re- 
mains insufficient to reconcile the expected signal with 
the observed cross-correlations. 

In this paper, we present the detection of cos- 
mic magnification using a large, uniform sample of 
photometrically-selected SDSS galaxies and quasars. 
Contrary to previous results, we find an excess of bright 
(g < 19) quasars around galaxies and a deficit of fainter 
quasars, matching the expected variation with quasar 
number count slope. In addition, the amplitude of the 
signal and its angular dependence is, for the first time, in 
agreement with theoretical predictions. Based on a num- 
ber of data quality and uniformity tests, we find that this 
detection is robust against possible sources of systematic 
error that may have plagued previous measurements and 
represents a genuine detection of magnification bias. The 
outline of the present paper is as follows: previews the 
basic weak lensing models for the expected signal and 
121 describes the galaxy and quasar data sets and cross- 
correlation estimators. 21 summarizes the results from 
magnitude-limited quasar samples as well as optimally- 
weighted measurements using the full quasar sample and 
compares them to the expected signals derived in £J2I Fi- 
nally, discusses the possible applications for further 
measurements using the SDSS and future large area sur- 
veys. 

2. MODELING MAGNIFICATION STATISTICS 

In this section we briefly describe the formalism of cos- 
mic magnification and introduce the notation that will be 
used below. Let no(/) df be the number of sources with 
a flux in the range [/, f + df) and n(/) df the correspond- 
ing number of lensed sources undergoing a magnification 
/i. We write the unlensed source counts as 



no(/)d/ = a r s(/) df. 



(1) 



where ao is some normalization factor and s(f) is the 
power-law slope as a function of flux /. The magnifica- 
tion effect will enlarge the sky solid angle, thus modifying 
the source density by a factor and at the same time 
increase their fluxes by a factor /i. These effects act as 
follows on the number of lensed sources: 

n(/)d/ = - n - — 
M V/V 



-2 if 



-«(///*) 



(1/ 



(2) 



If s does not vary appreciably over the interval [/, //z] , 
which is well satisfied if /i departs only weakly from unity, 
then 

n(/)d/ = M s(/) - 2 n (/)d/. (3) 

Expressing this as a function of magnitude, we recover 
the form appr opriate for a magnitude limited sample 
()Naravar]ll989]) : 

N(m) dm = ^ s(m) " 1 N (m) dm 

= ) u«M- 1 n ( m )dm, (4) 

The final form of the exponent (a(m) — 1) reflects the two 
distinct effects of magnification and how they interact to 
produce the signal observed on the sky: the amplification 
effect that varies as a function of the quasar magnitude 



and the dilution effect that is a constant regardless of 
magnitude. The combination of these two effects is the 
magnification bias. 

In the statistical context, magnification creates cor- 
relations between foreground and background popula- 
tions. In the weak lensing regime, i.e. if the conver- 
gence (k) and the shear (7) are small compared to unity, 
a first-order Taylor expansion of the magnification gives 
fi « 1+2 k. Using this approximation, a cross-correlation 
between magnification and foreground matter overdensi- 
ties can then be easily computed as a function of the mat- 
ter power spectrum. The formalism for mag nification by 
large- scale structures was first introduced bv lBartelmannl 
(1995) . Following the pres c riptio n and notation laid out 
in lJain. Scranton &: Shetbl l)2003[) . the expected magnifi- 
cation bias signal can be written as 



wgq(^i m) = 12ir 2 n M (a(m) - 1) X 



J d X dk kK.(k,e,x)Pg 



(5) 



= (a(m) — 1) x wq(6), 

where m is the magnitude of the sources, f^M is the cos- 
mological matter density relative to critical, \ is the co- 
moving distance, /C is the lensing kernel, and P gm (k, x) is 
the galaxy-dark matter cross-power spectrum. This for- 
mulation separates the expected signal into two pieces: 
wq{9), which contains all of the information about non- 
linear galaxy biasing and the redshift distributions; and 
{aim) — 1), which varies only with quasar magnitude 
and controls the sign of the expected signal. We note in 
passing that wq is closely related to the tangential shear 
induced by the galaxy-mass correlation, which is mea- 
sured in galaxy-galaxy lensing. In practice, one needs to 
consider quasars over a given magnitude range, in which 
a(m) may vary. In this case we have: 

w GQ (d) = (a-l}xw (e), (6) 



where 



1) 



JdmN(m)(a(m) — 1) 



(7) 



J dmN(m) 

The lensing kernel /C in Eq. EJis primarily a function of 
th e redshift distributions o f the g alaxies and quasars. As 
in I Jain. Scranton fc Shetbl (|2003f> . we model the redshift 
distributions as a combination of power-law and expo- 
nential cut-off: dN/dz ~ z a exp(— (z/z ) b ). The galaxy 
redshift distribution can be inferred from th e luminos- 
ity fu nctions measured by the CNOC2 survey l)Lin et alJ 
1999) with the appropriate apparent magnitude limits 
(see ®. For the quasars, photometric redshifts are com- 
puted for each quasar (see along with upper and 
lower redshift bounds and the probability that the quasar 
is within that redshift range. To model the redshift dis- 
tribution for each quasar sample, we assume a flat dis- 
tribution between the upper and lower redshift bounds 
for each quasar and weight according to the aforemen- 
tioned redshift probability. To first order the redshift 
distributions for all five quasar magnitude cuts are in- 
distinguishable, so we use the same fitted distribution to 
model each cut, making (a — 1) the only free parame- 
ter separating each magnitude bin. Fitting the redshift 
distributions to the chosen form, we find: 

^z 13 exp[-(z/0.26) 217 ] 

/ G 
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~z 256 exp [-(z/2.02) 12 - 76 ] , (8) 

Q 



where the redshift distribution for the quasars is limited 
to the range 1 < z < 2.2 using photometric redshifts (see 
SJSJ. Figure ^ shows the raw and fitted redshift distribu- 
tions. 

The only remaining piece of the theoretical calculation 
is the power spectrum. Since we are in the non-linear 
regime of gravitational collapse for the smallest angular 
bins and our foreground redshift distribution is broad, 
a thorough calculation of the expected signal would in- 
volve a matter-galaxy power spectrum calculated with an 
evolving halo occupation distribution (HOD). However, 
for the purposes of this paper, we are only interested in 
checking our measurements against a simple model of the 
expected signal, leaving the task of extracting the proper 
HOD behavior for future papers. With that in mind, we 
assume a simple HOD: 

(N) (A/) = 7V C + (M/M f , (9) 

where N c is unity for halo mass above 10 11 Ii~ 1 Mq and 
zero otherwise, Mo is 10 12 h~ 1 M Q and (3 is roughly unity. 
These parameters are approxi mately what one finds from 
semi-analytic galaxy codes (cf. iKravtsov et alJ2004|) and 
HOD fits to the SDSS s pectroscopic survey 2-point clus- 
tering measurements in lZehavi et alJ l)2004|) . 

We will use this formalism and a flat WMAP cos- 
m ology (fl M = 0.29 , fi A = 0.71, h = 0.72, n = 
1: iSpergel et al.ll200ffi to estimate u>o(9) and compare 
our measurements to the theoretical expectations. On 
small scales, non-linear magnification must be taken into 
account in order to obtain a more accurate modeling 
ijMenard et alJl2003t pFakada fc Hamanall2003|) . but this 
level of precision will suffice for our current case. A more 
detailed modeling which includes marginalization over 
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Fig. 1. — Galaxy redshift distribution from applying our 17 < 
r < 21 magnitude limit to the CNOC2 luminosity function and 
quasar redshift distribution inferred from quasar photometric red- 
shifts (solid lines). The fitted redshift distributions from Equa- 
tion [5] are shown with dashed lines. In all cases, the amplitude 
scaling is arbitrary. 



cosmological and redshift distribution parameters will be 
used in a future paper in order to constrain some of the 
model parameters. 

With this model in hand, we can test whether the 
measured signal is due to gravitational lensing in two 
ways: (i) we can test whether the amplitude of the cross- 
correlation properly varies as a function of magnitude, 
i.e. if WGQ(m) cx (a(m) — 1), where (a(m) — 1) is eval- 
uated according to Eq. and (ii) we can check if the 
angular variation of wgq(9) agrees with theoretical ex- 
pectations. Showing that the signal satisfies these two 
conditions is a robust test to demonstrate the lensing 
origin of the signal and a general lack of systematic con- 
tamination. 

3. DATA ANALYSIS 

The use of large, homogeneous s amples of galaxies 
and quasars ob served by the SDSS E mm "et~aLlll9M 
iFukugita et"aH 119961 iSmith et al] 12002ft improves on 
previous measurements of the cosmic magnification in 
several key ways. First , the SDSS provides accurate 
multi-color photometry llLupto n. Ounn fc Szalayl 1199 91 
llvezic et all l2004t iHogg et alJ 120011 iPier et alJ I2003T> 
over large areas of sky with tight control of the sys- 
tematic errors that could alter the observed density 
of sources on the sky (e.g., seeing variations, masks 
around bright stars, sky background subtraction prob- 
lems around bright galaxies). Moreover, the accuracy 
of the photometry is crucial, especially to determine the 
number of faint sources. If the required photometric ac- 
curacy is a few percent, CCD-based photometry is clearly 
superior to photographic plate data. For example, the 
2dF photometric accuracy is approximately 0.2 mag for 
objects with 17 < bj < 19.45. The corresponding incom- 
pleteness introduces an extra scatter for the density of 
sources on small scales and can thus mimic a signal. 

The second advantage of the SDSS comes from the 
multi-color photometry. As described below, consistent 
color-based selection over the full photometric survey al- 
lows us to select a larger (both in quasar numbers and in 
area), more uniform quasar sample than has ever been 
compiled for cross-correlation studies. This is critical 
both for minimizing Poisson errors on the measurement 
as well as avoiding systematic selection effects. In addi- 
tion, the multi-color photometry allows for the reliable 
estimation of photometric redshifts for quasars, remov- 
ing any redshift overlap between our galaxy and quasar 
populations. Given the fact that correlations due to in- 
trinsic clustering have a much larger amplitude than the 
ones expected from lensing, even a small fraction of back- 
ground sources at low redshift can give rise to a positive 
amplitude bias in the cross-correlation. 

3.1. The Data 

The data s et was drawn from th e third SDSS data 
release (DR3; lAbazaiian et al.1 12005). Before masking, 
this set covers roughly 5000 square degrees, the ma- 
jority of which is located around the North Galactic 
Cap. To limit our conta mination from systematic er- 
rors ijScranton et aljr2002|) in the photometric data, we 
imposed a seeing limit of l'/4 and an extinction limit of 
0.2 in the r band. We also included a mask blocking 
a 60" radius around bright galaxies (r < 16) and stars 
with saturated centers to avoid losing quasars due to lo- 
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cal fluctuations in sky bri ghtness and observing defects 
ijMandelbaum et all 120051) . The combination of these 
masks reduced our total area to ~ 3800 square degrees. 

With these systematic cuts, we can reliably perform 
star /galaxy separation using Bayesian methods to r = 21 
l)Scranton et alJfeOO^ . These selection criteria yielded 
13.5 million galaxies between 17<r<21ata den- 
sity of approximately one galaxy per square arcminute. 
For galaxies, we use countsjnodel magnitudes, while 
quasar magnitudes are given using psf counts magni- 
tudes ijStoughton et al.ll2002t these magnitudes are des- 
ignated as modelMag and psfMag in the SDSS database, 
respectively). In all cases, we de- redden the magnitudes 
to correct for Galactic dust extinction before applying 
the various magnitude cuts. For even modestly faint 
magnitudes (r > 18), the Petrosian magnitudes used in 
the SDSS spectroscopic sample ( and the previous SDSS 
galaxy-quasar measurements by iGaztanaeal 12003(1 can 
fluctuate with the local seeing, leading to a strong vari- 
ation (~ 25%) in apparent galaxy density with seeing 
for a magnitude-limited sample. The apertures used for 
countsjnodel magnitudes are convolved with the local 
PSF which makes them much more robust against seeing 
variations. For the seeing range between 0'.'85 and 1'.'4, 
the observed galaxy density as a function of local seeing 
is constant for our magnitude cut. As mentioned in §21 
applying these apparent magnitude cuts to the CNOC2 
luminosity functions yields a mean redshift for this mag- 
nitude limited sample of z ~ 0.3, with the maximum 
redshift of the sample near z ~ 0.75 

The quasar data set was generated usin g kernel density 
estima tion (KDE) methods described in IRichards et all 
1(2004(1 . Although our quasar sample is drawn from 
DR3, the selection method is identical to the one 
IRichards et alJ 1(2004(1 applied to the DR1 data set. The 
KDE method is a sophisticated extension of the tradi- 
tional color selection technique for identifying quasars. 
In this implementation, two training sets, one for stars 
and one for quasars are defined. Then the colors for each 
new object are compared to those of each object in each 
training set and a 4D Euclidean distance is computed 
with respect to the objects in each training set. New ob- 
jects are then classified in a binary manner (quasar/star) 
according to which training set has a larger probability of 
membership. This technique allowed clean separation of 
relatively low redshift (z < 2.5) quasars from the stellar 
locus in 4 dimensional color space, producing a catalog of 
225,000 quasars down to a limiting magnitude of 21 in the 
g band with greater efficiency and comple teness than the 
SPSS spectroscopic targe ting algorithm (Ric hards et alJ 
l2002tlBlanton et al.l2003|) . After masking, the total pop- 
ulation was reduced to 195,000 quasars. 

In addition to find ing quasars, we a pplied photometric 
redshift techniques fWeinstein et alJ[2004() to filter out 
low redshift quasars which might be physically associated 
with our foreground sample. Given the broader features 
and larger redshift range for quasars relative to those of 
galaxies, the photometric redshift errors for the quasars 
are generally somewhat asymmetric. Rather than esti- 
mate a Gaussian redshift error, we used an upper and 
lower redshift bound along with the likelihood that the 
redshift was within those bounds. To prevent redshift 
overlap with the galaxies, we required that the upper 
and lower bounds were within the range 1 < z < 2.2 and 



weighted each quasar according to the aforementioned 
redshift likelihood for both the number count and cross- 
correlation measurements. 

3.2. Measurement 

The expected lensing signal for magnification bias is 
generally dominated on small scales (< 0.01°) by Pois- 
son noise and falls below the noise on scales larger than 
1°. To cover this full range (and beyond), we used two 
estimators. For angular scales below 0.1°, we used a 
pair-ba sed estimator similar to the Landy-Szalay esti- 
mator l(Landv fc Szalavjll993() . but modified for a cross- 
correlation: 

m (GQ) ~ (RgQ) - (GRq) + (RgRq) nn , 
Wgq{9) = (RM ' (10) 

where (GQ) is the number of galaxy-quasar pairs sepa- 
rated by angle 9, (RgRq) is the number of pairs of ran- 
domized galaxy and quasar positions separated by 9, etc. 
To limit the Poisson noise in our estimation of (RgQ), 
(GRq) and (RgRq), we generated 50 random points for 
each galaxy and quasar. 

As we move from small to large scales, the estimator in 
Equation^Jbecomes progressively less and less efficient; 
As the angular scale increases, a progressively larger area 
much be searched for suitable pairs, and the estimator 
in Equation 1101 becomes less and less efficient, increas- 
ing computation time. Thus, for anglar bins larger than 
0.05°, we used a pixel-based estimator. Calculating the 
fractional galaxy and quasar over-densities (Sg and 6q, 
respectively), wqq is given by 



wgq(9) 



fifi 



(11) 



where we sum over all pairs of pixels separated by an- 
gle 9 and fi is the fraction of pixel i that remains after 
masking. This angular split roughly divides the total 
computation time for all of the various sub-samples (see 
|1} equally between the large and small angle estimator 
codes. 

For both estimator s, we used 30 jack-knife samples 
l(Scranton et al]l2002fl to generate errors, allowing us to 
combine the large and small angular measurements (in- 
cluding the single overlapping angular bin) to generate a 
coherent covariance matrix (C(9 a ,9p)), 



C(6 a ,6p) — 



N 



N - 1 



Eill ( w GQ,i(9a) 



W G Q(9 a )) X 

wgq(9p)), 



(12) 



where N is the number of jack-knife samples and 
wgq(9o) is the average value of u>GQ,i(#) for all N sam- 
ples. Equation 1121 measures the variance directly on 
the sky, so it should capture the contribution from the 
cross-correlation as well as the galaxy and quasar auto- 
correlations. We expect the errors to be dominated by 
Poisson noise, with subdominant terms coming from cos- 
mic variance as well as lensing by foreground structure 
not contained in our galaxy sample. To increase our sen- 
sitivity at small angles where the signal is most interest- 
ing, we employed a hybrid logarithmic binning scheme. 
For the angular decade running from 0.001° to 0.01°, we 
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Magnitude (g) 

Fig. 2. — Weighted number counts per magnitude per square de- 
gree in the 5 magnitude bins for quasars with photometric rcdshift 
between 1 < z < 2.2. 



used 3 logarithmically spaced bins, 4 bins for 0.01° to 
0.1°, 5 bins for 0.1° to 1°, etc. This improved the signal- 
to-noise on small scales at the expense of generating a 
slightly larger off-diagonal elements in the covariance ma- 
trix than that produced by a straight logarithmic binning 
system. However, the covariance matrices remained in- 
vertible in all cases with no degenerate modes, allowing 
us to use them for significance testing and curve fitting 
with no complications. 

4. RESULTS 

4.1. Lensing Origin of the Signal 

In order to investigate the lensing origin of the sig- 
nal, we first measured the quasar-galaxy correlations 
wgq(^, m) as a function of quasar magnitude in a given 
band. In Figure[21we show the number counts of quasars 
as a function of magnitude in the g band. We separated 
the g-selected quasar sample into five magnitude ranges 
and estimated the corresponding value of (a — 1} from 
power law fits in each bin. The results are presented in 
Table l4~Tl As can be seen, the values of (ct—1) are greater 
than zero for the three brighter magnitude bins. There- 
fore, due to the magnification bias, we expect to find 
an excess of such quasars in the vicinity of foreground 
lenses. In a similar way, we expect a deficit of quasars 
with g > 20. 

The corresponding quasar-galaxy correlation functions 
in each magnitude bin are shown in Figure [3] As 
expected, the brightest quasar sample with the steep- 
est number count slope showed the strongest positive 
cross-correlation, with the signal amplitude dropping and 
eventually changing to an anti-correlation as (a — 1) de- 
creases. At large angles, we see a flat signal consistent 
with zero for all five quasar magnitude bins. This first 
test verified qualitatively that the measured signal satis- 
fied the first of the two criteria described at the end of [j2I 
amplitude variation as a function of quasar magnitude. 

We can now quantify this agreement by using the 



Magnitude 


(a-1) 


17 < g < 19 


0.95 


19 < g < 19.5 


0.41 


19.5 < g < 20 


0.07 


20 < g < 20.5 


-0.24 


20.5 < g < 21 


-0.50 



TABLE 1 
The weighted mean value of 
a(m) — 1 = 2.5dlogNo(m)/dm — 1 OBTAINED FROM POWER LAW 
FITS IN DIFFERENT MAGNITUDE BINS IN THE g BAND. THE VALUES 
OBTAINED IN ALL FIVE BANDS ARE PRESENTED IN TABLE |^] 



-0.02 
0.O4 



-0.02 
0.04 



17 < g < 19 -- 




19.5 < g < 20 



ITU 



20.5 < g < 21 - 



1 i i 1 i 



19 < g < 19.5 - 



4- 



20 < g < 20.5 



"Ft 



O.Ol O.l 

6 (degrees) 



0.01 0.1 

6 (degrees) 



Fig. 3. — Measurements of wqq(9) as a function of quasar g 
band magnitude. Error ba rs a re the la errors based on the jack- 
knife covariance (Equation [HJ. The dark solid curve is the fitting 
model and the light dashed curve is the expected curve from the 
number counts {a — I). The shaded region indicates the la range 
on the fitted value of {a — I). Fitted and expected values for each 
magnitude bin are given in Table For angular scales larger than 
1 degree, the measurements were consistent with zero in all five 
magnitude bins. 



model given in arid the covariance matrices measured 
using Eauationll2l (see Figure^} to fit the measured data 
points and estimate the value of (a — 1} in each magni- 
tude bin. These fits are shown with the solid black line in 
Figure|3]and the one-sigma uncertainty by the shaded re- 
gion. The value of the parameter {a — 1} obtained in this 
manner can be compared to the one directly measured 
from the quasar number counts. The expected measure- 
ment based on the quasar number counts is shown by the 
dashed red line. For all five ^-selected magnitude bins, 
we find agreement between the fitted and measured val- 
ues of {a — 1} as a function of quasar magnitude. This 
demonstrates that the behavior of the signal quantita- 
tively follows both the amplitude and the angular varia- 
tions expected from magnification bias. 

We repeated similar measurements using magnitude- 
limited samples in each of the other four SDSS bands. 
The quasar number counts in all five bands are given in 
Figure [S] As mentioned above, the quasars were magni- 
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6 (degrees) 



Fig. 4. — Normalized covariance matrix for the 17 < 
g < 19 measurement of wqq(9). For a given pair of angu- 
lar bins (i and j), the normalized covariance matrix is given by 
C(6i, 0j)/ y/C(6i, e z )C(6j , 6j). The level of correlation between 
angular bins is roughly consistent for all of the measurements of 
w>gq(0). 
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Fig. 5. — Number counts per magnitude per square degree in 
the 5 filters for quasars with photometric redshift between 1 < 
z < 2.2. The original sample is magnitude limited in the g filter. 
The effective color cuts resulting from the quasar selection lead to 
incompleteness in the other filters at the faint end. 



tudc limited in g. For the other bands, the combination 
of effective color cuts for the sample and intrinsic scat- 
ter led to strong incompleteness for magnitudes fainter 
than 20. For the purpose of separating the sample into 
magnitude bins, the turn-over point set the faintest limit 
for each band. The results for these measurements are 
summarized in Table As with the (/-selected measure- 
ments, cross-correlations in the other four filters found 
qualitative and quantitative agreement with the expected 
magnitude bias variation with (a — 1). 
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Fig. 6. — a — 1 and (a — l) 2 as a function of magnitude in the 5 
filters for quasars with photometric redshift between 1 < z < 2.2. 
For the u, r, i, and z bands, incompleteness in the number counts 
at the faint end causes a — 1 to rapidly diverge, making a direct 
measurement of the expected anti-correlation in these filters very 
difficult. 



4.2. Optimal Stacking and Detection Significance 

Having shown that the signal follows the theoretical ex- 
pectations as a function of magnitude, we combined these 
measurements to quantify the significance of the global 
detection of cosmic magnification. Instead of separating 
the signal into five magnitude bins, we measured the sig- 
nal integrated over all magnitudes weighted with differ- 
ent powers of (a(m) — l) n . In Figure|Sl we show the num- 
ber count relations in each of the five SDSS filters and 
Figure plots the corresponding values of (a(m) — l) n 
for n — 1 and 2. These plots are made by measuring 
{a — 1} in narrow magnitude bins over the full range and 
then interpolating over the bins with a cubic spline. 

• Mean correlation function: by simply averaging the 
signal from all quasars, i.e. considering the case n = 1, 
we recover Equation^! where (a — 1) is given by integra- 
tion over the full magnitude range of the sample. The 
effects of bright and faint quasars generally canceled each 
other, resulting in the small values of (a — 1} found in 
Tabled Figure □ shows the results for all five SDSS fil- 
ters, along with expected and fitted curves for {a — 1). 
As with the magnitude-selected samples, we see gener- 
ally good agreement between the expected and observed 
signals. Note that all the data points are below the one- 
percent level. 

• Optimal correlat i on fu nction: as shown by 
iMenard fc Bartelmannl l)2002j) . using n = 2 (i.e. look- 
ing at the second-order moment of the signal as a func- 
tion of magnitude) optimally weights the expected lens- 
ing signal. This maximizes the S/N of the detection since 
the signal is weighted proportionally to the expectations. 
With the extra factor of a(m) — 1, the expected signal 
is: 
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Fig. 7. — Same as Figure |3 but over the full magnitude range 
in each wavelength band. The range on the y-axis is the same as 
in Figure |21 to provide an easier comparison of the relative ampli- 
tude for samples where we do and do not expect a strong lensing 
signal. Fitted and expected values for (a — 1} are given in Ta- 
ble 1^1 The correlation between angular bins for each measurement 
is consistent with that shown in Figure |1] to first order. 



o.oi o.i l 

9 (degrees) 

Fig. 8. — Same as Figure 171 but for iUqq™"' (9) , which uses 
quadratic weighting to enhance the lensing S/N (Equation 1131 . 
Fitted and expected values for (a — 1)e are given in Table 121 The 
correlation between angular bins for each measurement is consis- 
tent with that shown in Figure|I]to first order. 



JdmNjm) (a(m)-l) 2 

= r - ; 777 — n x Wo(") 

J dmN(m) v ; 

= (a-l) E xw Q (9). (13) 

The corresponding signal can be measured by weight- 
ing each quasar by a factor of a(m) — 1 and by calcu- 
lating the cross-correlation in the manner described by 
EauationsllOlandllll Rather than largely counter-acting 
each other as seen in Figured the positive and negative 
correlations from the bright and faint end of the quasar 
number counts now act in concert and benefit from the 
statistical power of the entire quasar population. The 
corresponding results are presented in Figure |S] Once 
again, we find a very good agreement between the model 
and the observations for all five bands. 

Using this optimally- weighted correlation function and 
the associated covariance matrix, we can assess the signif- 
icance of our detection. By comparing the corresponding 
X 2 values of Wqq'™^) against the null for 18 angular 
bins, we detect the signal at 4.1a, 8.1er, 4.8cr, 5.4<r and 
4.8c in the it, g, r, i and z bands respectively. If we con- 
sider only the angular scales < 1°, as shown in Figure 
the significance of the detection remains nearly the same. 

Given the high S/N provided by the optimally- 
weighted estimator, we can compare the angular vari- 
ation of the measured signal to theoretical expectations. 
As can be seen in FigJHI we find consistency from 0.3' 
to 1 deg, i.e. over more than two orders of magnitude 
in scale. This allows us to validate the second criterion 
from 321 the match of the angular variation of the signal 
with the predicted cross-correlation function. Consider- 
ing an effective redshift of z = 0.3 for the foreground 
galaxy population, we find that the detected magnifica- 
tion signal probes scales ranging from ~ 60/i _1 kpc to 10 
Mpc. 



Observing both the first and second moments of the 
lensing signal as a function of magnitude to give the ex- 
pected behavior as a function of the observed values of 
a(m) — 1 is an excellent indication that we are observing 
the signal originating from gravitational lensing. 

4.3. Systematics 

As seen in ^ accurate measurements of galaxy and 
quasar number counts can suffer from a number of biases: 
seeing variations, stellar contamination, dust extinction, 
redshift overlap, etc. To verify that our measurements 
are not affected by these effects, we have performed a 
number of checks. 

To test for stellar contamination in our sample of 
quasars, we cross-correlated stars in the 17 < r < 21 
range with the g band selected quasars, wsq(9), both 
with and without optimal weighting. Unlike galaxies or 
quasars, the local stellar density is not well approximated 
by the global mean density. As a result, we do not ex- 
pect (and do not observe) a null correlation between stars 
and quasars (or stars and galaxies). Rather, our ob- 
served wsq(9) was consistent with the observed galaxy- 
star cross-correlation, both of which are consistent with 
a very small (~ 1%) level of stellar contamination. More 
importantly, when we optimally weighted u>sq(#), the 
signal was consistent with zero at all angular scales, as 
would be expected for a stellar density independent of 
a(m) — 1. This was in marked contrast to Wqq 1 ™'^), 
detected at 8cr. Cross-correlations with local seeing pro- 
duced similar results. 

We also tested the robustness of the photometric red- 
shift likelihood by applying a series of redshift likelihood 
thresholds (i.e. requiring that the probability that the 
quasar was within the upper and lower redshift ranges 
specified by the quasar photometric redshift algorithm 
was above a given value: 50%, 60%, 70%, etc.). In all 
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cases, the resulting galaxy-quasar cross-correlation was 
consistent with measurements made with no threshold, 
verifying that our signal was not dominated by low prob- 
ability outliers. Finally, a cross-correlation with low red- 
shift quasars produced a large amplitude, positive sig- 
nal. However, this last point was sensitive to a restrictive 
cut on the quasar redshift probability (> 80%) due to a 
strong shift in the probability distributions for quasars 
below z ~ 1; higher redshift quasars tended to have much 
higher redshift probabilities (~ 0.8) than low redshift 
quasars (peaks around 0.5 and 0.8). 

Next, we checked for possible cont amination by large, 
bright galaxies. As described in Ma ndelbaum et all 
( 2005), the estimation of the density of faint sources 
around bright extended objects can be biased induced 
by uncertainties in the sky subtraction. We do not ex- 
pect the sky subtraction issues to be as significant since 
our quasars are point sources, but, as described in fJU we 
applied a 60" mask around all r < 16 galaxies. Measure- 
ments with and without these masks were identical, but 
we included the masks in our final analysis to avoid any 
unforeseen effects. 

Finally, a bias that is not related to the data analysis 
but that might be intrinsically present is extinction by 
dust. Indeed, the presence of dust around galaxies is ex- 
pected to redden and extinct background sources. So far, 
the amount of dust on large scales has been poorly con- 
strained and its effects on measurements of quasar-galaxy 
correlations has been uncertain. However, our analysis 
indicates that the deficit of quasars due to dust extinction 
is subdominant to the density changes induced by gravi- 
tational lensing. Indeed, the fact that the measured sig- 
nal for the first and second moments behaves as expected 
as a function of the slope number counts, a(m) — 1, indi- 
cates that the signal might not be contaminated by other 
sources than gravitational lensing. Biases like dust ex- 
tinction or the above-mentioned effects are not expected 
to scale proportionally to a(m) — 1 and would there- 
fore affect the first and second moments of the signal 
in different ways. This would prevent the simultaneous 
agreements found above. Therefore we conclude that our 
current measurements are not significantly affected by 
biases, but a parallel effort is underway to quantify the 
reddening effects of the lensing galaxies more precisely. 

5. DISCUSSION & FUTURE APPLICATIONS 

In this paper we have presented a detection of cos- 
mic magnification obtained by cross-correlating distant 
quasars and foreground galaxies. Using data from ap- 
proximately 3800 square degrees of the SDSS photo- 
metric sample, we have cross-correlated the position of 
~ 200,000 photometrically selected quasars and large- 
scale structures traced by over 13 million galaxies, and 
we have detected a signal on angular scales from 20" to 
1 degree at high significance. 

The magnification bias due to weak lensing gives rise 
to an excess or a deficit of background sources in the 
vicinity of foreground galaxies depending on the value 
of the power-law slope of the source number counts: 
a(m) = 2.5dlogNo(™)/dm. Our measurements of the 
galaxy-quasar cross-correlation function exhibit the ex- 
pected behavior: bright quasars, with steep number 
counts, appear to be in excess around galaxies and large- 
scale structures, and faint quasars with shallow number 



counts are seen to be in deficit. On all scales, we find 
WGQ( m ) tx (a(m) — 1) in the five SDSS bands, as ex- 
pected. 

We have measured the first and second moments of 
the signal as a function of quasar magnitude and the re- 
sults are in very good agreement with what is expected 
from the magnification bias: depending on the band, the 
first moment gives an amplitude consistent with zero or 
smaller than ~ 5 x 10 -3 , as a result of the opposite ef- 
fects arising from the bright and the faint quasars. The 
second moment, which turns out to be the optimal signal 
estimator, exhibits a strong signal detected at > Aa in all 
five SDSS filters and reaching up to 8.1a in the g-band. 
The quasars are magnitude selected in the g-band, giv- 
ing us the largest sample in this band (other bands lose 
quasars at the faint end due to the effective color cuts in 
these bands), so the difference in the signal-to- noise in 
the other bands relative to g is unsurprising. Using this 
estimator we find the angular dependence of the signal 
to be in very good agreement with theoretical estima- 
tions of lensing by large-scale structures. Our measure- 
ments probe physical scales ranging from 60ft.~ 1 kpc to 
10/i _1 Mpc at the mean lens redshift. Since we do not 
expect the biases from systematic errors to scale pro- 
portionally to a(m) — 1, the simultaneous agreement of 
the first and second moments of the signal as a function 
of magnitude indicates that these systematic biases (in- 
cluding dust extinction) do not significantly affect our 
measurements . 

The SDSS quasar and galaxy samples used in our anal- 
ysis are significantly larger, more uniform and better 
characterised than any data sets used for this measure- 
ment previously. We have shown that biases including 
seeing variations, stellar contamination, sky subtraction 
issues and errors in the photometric redshifts are well 
controled and do not significantly affect the measure- 
ments. Whereas previously claimed detections reported a 
signal much larger than theoretical predictions, our mea- 
surement shows, for the first time, the expected ampli- 
tude and angular dependence for the standard cosmolog- 
ical model and a realistic galaxy biasing. As such, we 
conclude that the disagreement between theoretical pre- 
dictions and previous measurements was most likely due 
to larger systematic effects in these data sets which could 
not be adequately controled. 

The successful detection of cosmic magnification opens 
the door to a number of applications: as mentioned in 
cosmic magnification is a function of the first moment 
of the galaxy halo occupation distribution (HOD) on all 
angular scales. Conversely, the galaxy auto-correlation 
function, w(6), is a strong function of the second mo- 
ment of the HOD. Thus, measuring these two quantities 
for the same sample of galaxies will provide us with con- 
strains on both moments, and therefore probe the scales 
on which the galaxy biasing becomes stochastic. Such an 
analysis can then be carried out as a function of galaxy 
type, redshift, etc. and provide interesting constraints on 
our understanding of galaxies and large-scale structures. 

As noted in ^ our measurements of cosmic mag- 
nification constrain the projected galaxy-mass correla- 
tion in much the same way as galaxy-galaxy lensing, 
although that method is based on galaxy shapes and 
measurements of shear. This complementarity is a par- 
ticularly useful cross-check since the dominant sources 
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of systematic error for the two methods are different 
(PSF anisotropy for the galaxy-galaxy shear measure- 
ments vs. photometric calibration for the magnification 
bias). Furthermore, using quasars as sources, cosmic 
magnification allows for probing lensing at higher red- 
shifts: the sources used for SDSS galaxy-galaxy lensing 
studies are used as lenses for measurements of quasar- 
galaxy correlations. Finally, as is the case for cosmic 
shear, higher-order statistics can also be investigated in 
the context of l ensing-induced quasar-galaxy correlations 
ijMenard et al.ll2C)03jV . We also note that the techniques 
used for efficient quasar selection are readily applicable 
to next generation of large, multi-band surveys. Cos- 
mic magnification is therefore an excellent complement 
to planned cosmic shear surveys. 

In a future work we will use measurements of quasar- 
galaxy correlations to generate the first constraints on 
the galaxy HOD and cosmological parameters obtained 
from cosmic magnification. Likewise, projects are under- 
way to use cosmic magnification to measure the extent of 
galaxy dust halos as well dark matter halo ellipticities. 
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Scranton et al. 



TABLE 2 

Quasar counts, (a — 1) from the number count slopes in 

EACH MAGNITUDE BIN, (a - 1) FROM FITTING OBSERVED Wgq(O) 
TO THE WEAK-LENSING MODEL GIVEN IN EQUATION [H] AND 
DETECTION SIGNIFICANCE FOR THE OPTIMAL ESTIMATOR IN 
EACH FILTER. FOR THE OPTIMAL ESTIMATOR, THE FULL 

MAGNITUDE RANGE IS USED AND THE VALUES OF (a — 1} ARE FOR 

(a — 1)e as given in Equation El Detection significance 

IS GIVEN IN MULTIPLES OF STANDARD DEVIATION (cr). 



Magnitude Limit 


Quasar Counts 


(a - 1) 


Fitted (a — 1) 


a 


17 < u < 19 


6774 


+0.97 


+ 1.63 + 0.70 




19 < u < 19.5 


9001 


+0.42 


+0.19 + 0.95 




19.5 < u < 20 


16648 


+0.07 


-0.27 + 0.39 




20 < u < 20.3 


14824 


-0.11 


+0.04 + 0.55 




20.3 < u < 20.6 


18524 


-0.37 


-0.08 + 0.34 




17 < u < 20.6 


65610 


+0.05 


+0.09 + 0.19 








-1-0 1 8 


in i q -|- n DQ 


4.1 


17 < g < 19 


8054 


+0.95 


+ 1.53 + 0.57 




19 < g < 19.5 


10312 


+0.41 


+0.49 + 0.81 




19.5 < g < 20 


18148 


+0.07 


-0.06 + 0.33 




20 < g < 20.5 


28751 


-0.24 


-0.12 + 0.36 




20.5 < g < 21 


39567 


-0.50 


-0.32 + 0.19 




17 < g < 21 


104683 


-0.12 


-0.02 + 0.14 




Optimal 


104683 


+0.22 


+0 20 ± 05 


8.1 


16 < r < 18.5 


4212 


+1.19 


+ 1.84 + 0.79 




18.5 < r < 19 


6120 


+0.65 


+ 1.36 + 0.86 




19 < r < 19.5 


12101 


+0.35 


-0.30 + 0.60 




19.5 < r < 20 


21141 


+0.03 


+0.14 + 0.25 




20 < r < 20.3 


18137 


-0.25 


-0.11 + 0.34 




16 < r < 20.3 


61596 


+0.16 


+0.37 + 0.17 




Optimal 


61596 


+0.20 


+0 18 ± 12 


4.8 


16 < i < 18.5 


5609 


+1.08 


+ 1.66 + 0.81 




18.5 < i < 19 


7813 


+0.56 


+ 1.06 + 0.96 




19 < i < 19.5 


15236 


+0.23 


-0.10 + 0.76 




19.5 < i < 20 


26173 


-0.06 


+0.21 + 0.33 




20 < i < 20.2 


17687 


-0.51 


-0.11 + 0.57 




16 < i < 20.2 


72391 


+0.05 


+0.26 + 0.14 




Optimal 


72391 


+0.21 


+0.24 + 0.13 


5.3 


16 < z < 18.5 


5812 


+1.00 


+ 1.58 + 0.77 




18.5 < z < 19 


8047 


+0.62 


+0.97 + 0.85 




19 < z < 19.5 


16056 


+0.29 


-0.05 + 0.37 




19.5 < z < 19.8 


15240 


+0.13 


-0.02 + 0.39 




19.8 < z < 20.1 


20177 


-0.27 


-0.25 + 0.49 




16 < z < 20.1 


65207 


+0.19 


+0.31 + 0.18 




Optimal 


65207 


+0.21 


+0.25 + 0.11 


4.8 



